NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


1 9960223  077 

THESIS 


CALIBRATION  OF  A  COHERENT  ACOUSTIC 
SEDIMENT  PROFILER  (CASP) 

AND  ANALYSIS  OF  SEDIMENT 
DISTRIBUTION  FROM  DUCK94 

by 

Karen  Marie  Kohanowich 
September  1995 


Thesis  Advisor:  Timothy  P.  Stanton 

Co-Thesis  Advisor:  Edward  B.  Thornton 

Approved  for  public  release;  distribution  is  unlimited. 

jmc  QtjALrrr  mapBcrsD  8 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  239-18  298-102 


Approved  for  public  release;  distribution  is  unlimited. 


CALIBRATION  OF  A  COHERENT  ACOUSTIC  SEDIMENT  PROFILER  (CASP)  AND 
ANALYSIS  OF  SEDIMENT  DISTRIBUTION  FROM  DUCK94 


Karen  M.  Kohanowich 
Lieutenant  Commander,  United  States  Navy 
B.S.,  Vanderbilt  University,  1982 


Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 


MASTERS  OF  SCIENCE  IN  METEOROLOGY 
AND  PHYSICAL  OCEANOGRAPHY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1995 


Author: 


Approved  by: 


^  p 

Timothy  P.  Stanton,  Thesis  Advisor 


Edward  B.  Thornton,  Second  Reader 

Robert  H.  Bourke,  Chairman 
Department  of  Oceanography 


111 


ABSTRACT 


The  Coherent  Acoustic  Sediment  Profiler  (CASP)  uses  four  acoustic  beams  at  two 
frequencies  (1.3  and  5.3  MHz)  to  measure  the  backscatter  from  sediments  and  infer 
concentrations  in  1.68  cm  bins  to  one  meter  range.  It  also  measures  three  velocity 
components  and  sediment  concentration  in  the  same  intersecting  volume  (0(1  cm^))  at 
a  sampling  frequency  of  36  Hz.  Empirical  calibration  of  the  sediment  profiling  capability 
of  the  CASP  has  been  performed  with  sediment  sampled  from  the  nearshore  zone  at 
Duck,  N.C.  The  sediment  is  characterized  by  a  lognormal  distribution  with  a  mean  grain 
size  diameter  of  0.18  mm,  and  standard  deviation  of  1.2.  The  1.3  MHz  sonar  measures 
sediment  concentrations  of  0(1x10"^  kgm'^)  to  25  kgm and  the  5.3  MHz  sonar  measures 
concentrations  of  0(1x10'^  kgm'^)  to  8  kgm■^  both  with  an  accuracy  of  10%.  A  grain 
size  discrimination  capability,  based  on  comparison  of  measurements  by  the  two 
frequencies  and  the  assumption  of  a  lognormal  sand  size  distribution,  can  detect  variations 
of  ±  0.05  mm  mean  geometric  radius  between  the  calibration  sand  and  the  sediment 
measured  in  the  field.  Field  data  from  a  day  characterized  by  narrow  banded,  nonlinear 
wave  forcing  during  the  DUCK94  experiment  show  that  episodes  of  high  sediment 
concentration  correlate  well  with  the  maximum  onshore  velocites  within  2  cm  of  the 
bottom,  but  decorrelate  above  this  level.  Data  is  ensemble  averaged  with  a  phase 
mapping  technique  to  illustrate  this  correlation.  Sediment  flux  calculations  reveal  strong 
shoreward  transport  at  the  boundary  layer,  and  weaker  seaward  transport  higher  in  the 
water  column. 
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I.  INTRODUCTION 


The  ability  to  model  sediment  transport  in  the  surf  zone  is  important  to  both  scientists 
and  engineers.  For  the  military,  it  is  an  integral  concern  for  the  planning  of  amphibious 
operations.  The  mechanisms  responsible  for  this  sand  transport  under  breaking  waves, 
however,  are  still  poorly  understood.  Among  the  most  elusive  problems  is  a  satisfactory 
description  of  the  small  scale  processes  in  the  boundary  layer,  a  deficiency  in  large  part  due 
to  a  lack  of  sufficient  spatial  and  temporal  resolution  of  sediment  flux  measurements. 

A  variety  of  instruments  have  been  used  to  measure  sediment  concentrations  in  the 
surf  zone,  including  physical  collection  of  sand  by  "sand  socks",  direct  siphoning,  and 
optical  and  acoustic  backscatter  sensors.  Physical  collection  techniques  are  limited  in 
volume  resolution,  extremely  limited  in  temporal  resolution,  and  can  disturb  the  in  situ  flow. 
Optical  backscatter  sensors  (OBS)  use  scattered  infrared  radiation  to  measure  suspended 
particle  concentrations,  and  are  often  coupled  with  nearby  current  meter  measurements  to 
estimate  flux.  Although  the  temporal  resolution  is  on  the  order  of  1 0  Hz,  the  sample  volume 
depends  on  the  concentration  of  the  sediments  being  measured,  and  as  a  result  is  poorly 
defined.  The  frame  on  which  the  instruments  are  mounted  in  the  field  also  tends  to  perturb 
the  flow,  and  the  dual  instrument  technique  cannot  measure  concentration  and  velocity  in  the 
same  volume.  Laser  Doppler  Velocity  (LDV)  measurements  have  been  limited  due  to  the 
high  cost  and  complexity  of  the  optical  systems,  and  are  usually  limited  to  two  velocity 
components.  The  sample  volume  of  the  LDV  is  only  100  pm  by  1  mm,  which  makes 
calibration  for  concentration  profiling  very  difficult.  Hot  film  anemometers  provide 
extremely  high  resolution  measurements  of  current  speed  within  the  bottom  boundary  layer, 
which  makes  them  ideal  for  the  study  of  turbulence  length  scales.  They  carmot  sense 
direction  of  the  flow  or  measure  sediment  concentrations,  however,  so  they  are  of  limited 
value  in  sediment  flux  studies.  Acoustic  instruments  are  able  to  make  non  invasive,  high 
resolution  measurements  of  sediment  concentration.  They  have  been  operated  successfully 
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in  the  field  for  several  years  (Hanes  and  Huntley  (1986),  Hay  (1985),  Thome  et  al.  (1991), 
Young  (1982))  and  typically  profile  sediment  concentrations  to  a  range  of  one  meter  with 
one  cm  and  one  second  resolution.  Sediment  concentrations  are  retrieved  from  backscatter 
signals  through  empirical  and/or  analytical  calibration  techniques.  Limitations  in  acoustic 
profiling  systems  occur  at  high  sediment  concentration  as  multiple  scattering  effects  begin 
to  occur.  Single  frequency  instmments  and  collocated  multiple  frequency  arrays  have  not 
provided  a  coincident  velocity  measurement. 

The  acoustic  instrument  used  in  this  research  was  designed  to  address  the  above 
limitations  (Stanton  1993).  The  Coherent  Acoustic  Sediment  Probe  (CASP)  consists  of  a 
forward  looking  1.3  MHz  beam  surrounded  by  three  5.3  MHz  beam  positioned  so  that  the 
four  beams  intersect  approximately  25  cm  in  front  of  the  instrument.  Using  coherent  doppler 
techmques,  three-dimensional  sediment  velocities  are  measured  in  this  intersection  volume. 
A  mean  grain  size  estimate  can  be  made  by  comparing  the  backscatter  characteristics  of  sand 
from  the  two  frequencies,  so  that  violations  of  homogeneous  sand  size  assumptions  can  be 
detected.  As  with  previous  instruments,  sand  concentration  is  also  measured  along  the  length 
of  the  beam.  During  the  DUCK94  field  experiment,  the  CASP  was  mounted  on  a  hydraulic 
arm  so  that  measurements  could  be  made  at  various  heights  above  the  bed. 

The  background  theory  for  use  of  acoustic  techniques  to  measure  nearshore  sediment 
profiles  is  presented  in  Chapter  II.  Chapter  III  describes  the  calibration  of  all  four  sonars 
for  measurement  of  sand  concentrations,  including  analysis  of  concentration  limits,  and  the 
capabilities  of  the  grain  size  discriminator.  Field  data  from  a  day  characterized  by  narrow 
band  wave  forcing  during  DUCK94  is  presented  to  demonstrate  the  capabilities  of  the 
CASP  in  Chapter  IV .  Utilizing  a  phase  mapping  technique,  vertical  sediment  distribution 
is  correlated  with  cross-shore  velocity  phase.  The  performance  of  the  grain  size 
discriminator  in  the  field  is  also  analyzed.  Chapter  V  presents  conclusions  and  areas  of 
future  research. 
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II.  BACKGROUND  THEORY 


Hay  (1991)  has  shown  that  the  pressure  amplitude  returned  from  an  acoustic  beam 
illuminating  sand  particles  is  directly  related  to  the  square  root  of  the  concentration  of  the 
particles.  Other  factors  involved  in  this  relation  are  the  decay  of  the  amplitude  of  the  signal 
due  to  spreading  of  the  sound  beam,  absorption  of  sound  energy  by  the  water,  and 
attenuation  of  the  signal  due  to  scattering  by  particles.  The  measured  acoustic  power  is  also 
affected  by  the  frequency  and  beam  characteristics  of  the  instrument,  and  by  the  composition, 
shape,  and  size  distribution  of  the  sand  particles.  There  are  essentially  two  approaches  to 
inverting  the  signal  return  to  provide  sand  concentration  estimates.  One  is  an  empirical 
method,  in  which  laboratory  derived  calibration  coefficients  describe  the  attenuation  and  the 
proportionality  between  the  square  of  the  acoustic  pressure  and  the  mass  concentration.  An 
assumption  necessary  for  this  method  is  that  the  size  distribution  of  the  sand  used  for  the 
laboratory  calibration  is  the  same  as  that  of  the  sand  measured  in  the  field.  The  second 
approach  is  analytical,  and  seeks  to  quantify  the  acoustic  response  in  terms  of  particle 
behavior.  This  method  uses  a  frequency  dependent  acoustic  form  factor,  which  accounts  for 
phase  shifts,  grain  diameter,  and  grain  size  distribution,  to  describe  the  total  scattering  cross 
section  and  scattering  attenuation.  The  absolute  backscatter  sensitivity  response  of  the 
system  is  then  scaled  to  a  given  lognormal  sediment  mass  concentration.  The  aim  of  the 
analytical  approach  has  been  to  exploit  frequency  dependent  acoustic  characteristics  in  order 
to  provide  estimates  of  sand  concentration  assuming  a  lognormal  distribution  with  an 
unknown  mean  grain  size. 

In  the  literature.  Young,  et  al  (1982),  and  Tamura  and  Hanes  (1986)  have  used  the 
empirical  approach.  Sheng  and  Hay  (1988),  Hay  (1991),  and  Hay  and  Sheng  (1992),  have 
described  a  series  of  analytical  models  for  multi-frequency  backscatter  and  the  response  to 
variable  mean  grain  size  and  lognormal  distribution  width.  Thome  et  al  (1991,  1993)  have 
used  a  combination  of  analytical  and  empirical  techniques. 
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Young  et  al  (1982)  first  calibrated  a  3  MHZ  acoustic  backscattering  profilometer  with 
three  sizes  of  sand  (120  gm,  200  gm,  and  450  gm)  in  a  laboratory  tank.  Concentrations  of 
from  .001  to  0.1  kg  m''’  were  compared  with  profiles  calculated  using  the  relation  ; 

(2.1) 

where  /,  is  the  backscattered  intensity,  or  voltage  squared,  K„  is  a  reference  concentration, 
corrects  for  spherical  spreading  with  range  r,  and  the  fixed  exponential  term  corrects  for 
attenuation  (attenuation  coefficient  =  1 .2  db/m).  Application  of  this  formula  to  backscatter 
profiles  confirmed  a  linear  relationship  for  each  of  the  sand  sizes  to  accuracies  within  a  factor 
of  two,  but  underestimated  the  concentration  of  the  largest  sand  particles. 

Tamura  and  Hanes  (1986)  conducted  further  calibrations  with  the  same  3  MHZ 
instrument.  They  used  full  range  and  sifted  samples  of  sand  collected  from  their  experiment 
site.  Stanhope  Lane,  in  concentrations  of  from  0.23  kg  m'^  to  5.7  kg  m'^.  They  determined 
the  semi-empirical  relationship  between  the  backscattered  intensity  /^(y)  and  the  sand 
concentration  C(y): 


0 

I,(y)=AC(y)*cxp(-J'(B^+B^C(y))dy)  (2.2) 

y 

where  the  constants  A,  a  function  of  the  acoustic  backscatter  cross  section,  B^ ,  the  loss  due 
to  water  absorption,  and  B^ ,  the  acoustic  forward  scattering  cross  section  of  the  sediment 
water  mixture,  were  determined  from  the  laboratory  experiments.  The  values  for  constants 
differed  between  the  full  range  and  sifted  sand  samples,  indicating  an  acoustic  response 
dependence  on  grain  size  and  distribution.  The  calibration  using  the  full  range  samples  was 
corrected  for  observed  attenuation  along  the  beam  and  was  accurate  to  within  ±10%. 
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Investigation  into  temporal  resolution  indicated  that  a  time  average  of  as  short  as  0.5  seconds 
still  produced  accurate  concentration  estimates.  The  analytical  method  strives  to  predict 
acoustic  backscatter  by  sediment  by  describing  all  aspects  of  sound  behavior  in  the  presence 
of  particles.  The  general  acoustic  equation,  which  describes  pressure  amplitude  of  sound 
waves  in  a  sample  volume,  is: 


(2.3) 

r 

where  P„  is  the  pressure  amplitude  at  the  reference  distance  ,  r  is  the  distance  from  the 
transducer  to  the  detected  volume,  and  is  the  attenuation  due  to  water  absorption.  To 
apply  this  relation  to  the  more  specific  problem  of  sand  concentration  measurement,  factors 
are  added  to  account  for  beam  spreading,  attenuation  due  to  scattering  by  particles,  loss 
through  the  measured  volume,  and  the  response  of  the  acoustic  beam  itself  to  the  particle 
population.  Through  such  applications,  Hay  (1991)  and  Hay  and  Sheng  (1992)  derived  the 
following  acoustic  model: 


V=S* 


—  *F[«(a)/„(x)]* 


N  P, 


£ 

sinh— 


-2  /  OLdr 


(2.4) 


where  V  is  the  output  voltage  of  the  receiver,  and  5  is  a  system  constant,  or  sensitivity, 
which  depends  on  the  receiver  directivity,  lobe  half  width,  reference  range  and  output 
pressure.  In  their  system,  the  time  variable  gain  (TVG)  of  the  receiver  compensates  for  the 
rg/r  spreading  loss  and  loss  due  to  absorption  by  water,  ^  the  mass  concentration  in 
kgm  \  and  is  the  particle  density.  F[/j(a)/„(x)]  is  an  empirical  function  which  includes 
the  effects  of  a  frequency  dependent  acoustic  form  factor,  fjx) ,  and  the  size  spectral  density, 
n(a).  Here,  x  is  the  nondimensional  particle  radius,  k^,  where  k,  is  the  wavenumber  of  the 


5 


sound  wave  in  the  fluid,  and  a  is  the  equivalent  radius  of  the  particle.  The  y/sinh^/^  term 
corrects  for  attenuation  across  the  detected  volume;  ^ -cx{a.^  +a^)  where  c  is  sound  speed, 
T  is  pulse  length,  and  a  is  attenuation  due  to  water  and  sediment,  respectively.  The 

r 

exponential  term,exp[-2y^dr]  describes  the  return  path  integrated  loss  due  to  scattering 
from  particles,  ° 

Of  the  above  variables,  the  sensitivity  S  can  be  measured  in  the  laboratory  for  each 
instrument.  The  form  factor  and  the  attenuation  due  to  scattering  by  particles  are  complex 
functions  of  the  compressibility  and  density  of  the  particles,  fluid  shape  of  the  particles,  and 
X.  The  form  factor  here  describes  the  acoustic  response  of  an  irregularly  shaped  sand  grain 
and  was  derived  semi-empirically  from  the  form  factor  for  spherical  quartz  grains.  This 
spherical  form  factor  has  been  described  by  Faran  (1962)  as  a  function  of  wavenumber, 
particle  radius,  phase  shift,  and  reflection  angle,  and  by  Gaunard  and  Uberall  (1983)  as  a 
function  of  wavenumber,  particle  radius,  the  velocity  of  sound  in  water,  compression  and 
shear  wave  velocities  of  the  sphere,  and  water  and  sphere  densities.  This  form  factor  for 
spheres  was  translated  into  one  for  irregular  particles  first  through  a  high  pass  model  for 
backscattered  intensity  developed  by  Johnson  (1977).  This  model  uses  a  polynomial  that  fits 
backscatter  behavior  in  the  Rayleigh  and  geometric  regions.  Sheng  and  Hay  (1988) 
expanded  the  high  pass  model,  based  on  Faran’s  (1962)  description,  by  including  angular 
dependence,  to  give: 


Kx^ 


(2.5  a) 


where  ^^=2/3  |  Yj^+YpCos6 


(2.5  b) 


6 


(K-K) 


.3(PrP.) 

C2p,+pj 


(2.5  c,d) 


where  k^,  ic„,  Ps,  pw  are  compressibility  and  density  of  sediment  and  water,  respectively. 
A  semi-empirical  fit  to  experimental  measurements  is  given  by  Sheng  (1991): 


toi= 


l+1.25x' 


Here,  \fj(x)^  \  is  the  theoretical  form  factor  for  a  rigid,  mobile  (finite  bulk  density)  spherical 
particle  with  the  density  of  quartz  (Hay  and  Mercer,  1985).  The  attenuation  coefficient  cCg 


is  described  as 


iVv 

a  = — 2j 
‘  2  ' 


where  N  is  the  total  number  of  particles,  and  2^  is  the  total  backscatter  cross  section  of  the 
particles,  which  describes  both  scattering  and  absorption  effects.  S*  can  be  represented  as 
an  explicit  function  of  the  spherical  form  factor  (Sheng  and  Hay  1988),  or  as  a  function  of 
X  and  compressibility  and  density  of  water  and  sediment  (Thome  et  al  1 993).  The  high  pass 
model  for  the  attenuation  coefficient  constmcted  by  Sheng  and  Hay  (1988)  is  expressed  as: 


a  ( - )= - 

M  [\+LKy+^'^-\ 


(2.8  a) 


(y^y;/3) 

where  K  =— — ^ — 
“  6 


(2.8  b) 
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and  I  is  an  adjustable  constant  5^1.  The  term  improves  the  fit  to  the  data  between 
Rayleigh  and  geometric  ranges.  From  the  known  values  and  from  laboratory  measured 
attenuation,  the  total  cross  section  can  be  calculated,  and  then  applied  to  the  forward 
problem  to  determine  attenuation  in  the  field. 

Hay  (1991)  used  three  frequencies  (1,  2.25,  and  5  MHZ)  to  examine  scattering 
behavior  in  the  laboratory.  To  avoid  nearfleld  complications,  the  instruments  were  aimed 
perpendicularly  into  a  jet  of  sediment,  and  corrections  were  applied  for  the  jet 
characteristics.  Results  for  [/I  |  calculations  show  very  good  agreement  with  theory  at  x  <1, 
but  give  values  of  L/L(x)  |  higher  than  expected  for  x  >  1.  The  curvature  of  the  fit  is 
somewhat  similar  to  theory  between  l<x<2,  but  departs  for  x>2.  Values  of  for  sand 
particles  was  calculated  using  direct  hydrophone  measurements  of  sound  propagating 
through  the  sediment  plume.  Although  subject  to  statistical  scatter,  the  data  clearly  fits  the 
modified  high  pass  model  better  than  the  rigid  sphere  theory.  Also  using  the  above  semi- 
empirical  expressions,  Thorne  et  al  (1993)  found  values  of  \fjx)  \  and  S,  calculated  from 
measurements  to  be  25%  more  than  those  predicted  by  theory. 

The  limits  of  the  linear  response  of  voltage  squared  to  mass  concentration  were  found 
to  be  24  kg  m’^  (1  MHz),  12  kg  m'^  (2.25  MHz),  and  8  kg  m'^  (5  MHz).  The  limit  for  the 
linear  response  of  attenuation  to  mass  concentration  was  determined  by  hydrophone 
measurements  to  be  24  kg  m'^  for  the  2.25  MHz  sonar  (and  its  4.5  MHz  harmonic),  and  for 
the  5  MHz  beam.  Attenuation  was  not  corrected  for  in  the  case  of  the  1  MHz  sonar. 
Measured  attenuation  values  were  used  to  correct  the  P  readings  for  the  2.25  and  5  MHz 
transducers.  This  linearized  the  2.25  MHz  values,  but  overcorrected  the  5  MHz  readings  for 
concentrations  greater  than  —10  kgm'^.  The  explanation  for  this  is  postulated  to  be 
averaging  of  a  non-linear  process;  i.e.  that  instantaneous  attenuation  must  be  applied  to  the 
instantaneous  backscatter  before  averaging  in  order  to  correctly  model  the  behavior. 
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Thorne  et  al  (1991)  used  3.0  and  5.65  MHz  sonars  in  a  combined  empirical  and 
analytical  approach.  The  transducers  were  aimed  longitudinally  through  a  cylinder  of  well 
mixed  sediment  in  the  lab,  and  corrections  were  made  for  nearfield  complications. 
Attenuation  was  calculated  following  Sheng  and  Hay  (1988),  but,  rather  than  calculate 
\fj(x)\,  a  calibration  constant  was  determined  from  a  single  range  and  concentration.  Overall 
accuracy  of  measurements  using  this  technique  was  ±20%.  They  point  out  that  attenuation 
along  a  beam  dominates  the  acoustic  response  to  increasing  concentrations  at  higher 
concentrations,  so  that  the  variation  in  <P^>  with  \Im  is  no  longer  linear.  In  fact,  the 
function  becomes  double  valued,  and  linear  empirical  calibrations  techniques  are  no  longer 
valid.  Their  investigation  also  demonstrated  that  a  10  second  average  of  10  Hz  signals 
clearly  presents  a  more  coherent  profile  than  the  signal  from  a  single  cycle.  Their 
calibrations  are  performed  using  the  10  second  average. 

Later,  Thome  et  al  (1993)  use  the  more  analytical  approach,  and  include  field  testing. 
They  illustrate  the  sensitivity  of  the  attenuation  integration  technique  by  performing  the 
integration  through  a  constant  sediment  concentration  of  2.2  kg  irf  .  A  5-10% 
underestimate  of  the  initial  concentration  led  to  an  order  of  magnitude  error  at  105  cm. 
Backscatter  results  from  field  deployment  compared  with  pumped  samples  from  the  same 
locations  yielded  agreement  to  -10%.  Errors  at  higher  concentrations  were  less  than 
laboratory  results,  presumably  because  the  concentration  was  not  as  consistently  high 
throughout  the  beam  as  it  was  in  the  laboratory. 

It  has  been  recognized  by  the  above  workers  that  the  use  of  more  than  one  frequency 
may  enable  grain  size  estimation  through  the  response  of  the  form  factor  and  thus  more 
accurate  concentration  inversion  procedures.  Hay  and  Sheng  (1992)  address  this  issue  using 
again  the  1,  2.25,  and  5  MHz  sonars.  When  the  assumption  is  made  that  scattering 
attenuation  is  small,  or  corrected  for.  Equation  (2.4)  can  be  expressed  as; 


F(X^o^)  /  =  1,2,5 


(2.9) 
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where  F(Xi  a  J  is  the  representation  of  F[n(a)X.  W]  under  the  assumption  of  a  lognormal 
sediment  distribution.  The  subscript  i  represents  the  different  frequencies,  Xi-k^g  (acoustic 
wavenumber  and  geometric  mean  radius)  and  =  the  standard  deviation  of  the  lognormal 
distribution.  The  ratios  of  signals  from  the  different  frequencies  can  then  be  written  to 
eliminate  the  dependence  on  M: 


F(X,op_^.r 


/^;,/=l,2,5y=l,2,5 


(2.10) 


Since  V  is  measured  and  S  is  calculated  independently,  the  left  hand  side  can  be  computed 
using  a  semi-empirical  expression  for  |/^  |  which  assumes  a  lognormal  sediment  (Hay  and 
Sheng  (1992)  Equation  (7)).  To  solve  for  both  unknowns,  and  o^,  two  independent  ratios, 
and  therefore  three  different  frequencies,  are  required.  Since  the  ratio  for  F^/F^  is  multi¬ 
valued  at  sand  diameters  less  than  100  pm,  and  the  slope  of  F^/F^  tends  to  0  over  100  pm, 
the  procedure  to  calculate  both  variables  is  often  indeterminate.  Therefore,  an  independent 
measurement  of  o^  was  made  on  sand  from  the  field  experiment  area,  and  used  for  the 
calculation.  Once  the  value  of  is  determined,  the  concentration  M  can  be  estimated 
from: 


p  Vl  Vf  i=\  for  a  <l(Wi\im 

- 1 - - ]  - - - 

2  slF\X^,a^)  sfF\X.,o^)  /=5 /or  100pm 


(2.11) 


Attenuation  is  addressed  with  a  bin  to  bin  integration  technique  that  incorporates  a  modified 
high-pass  model  for  introduced  by  Sheng  and  Hay  (1988).  Attenuation  for  the  1  MHz 
beam  is  assumed  to  be  negligible. 
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In  the  laboratory,  measurements  of  concentrations  up  to  23  kg  m‘^  of  particles  of 
from  90  to  430  pm  diameter  were  accurate  within  10%  using  the  three  frequency  algorithm. 
The  algorithm  underestimated  grain  size  by  10%  at  concentrations  above  0.8  kg  m'^,  and 
20%  at  lower  concentrations.  This  error  was  presumed  to  be  caused  by  weak  1  MHz  signal 
at  low  concentrations,  possible  deposition  of  coarser  materials  in  the  jet  mechanism,  and/or 
systematic  differences  between  the  measured  form  factor  for  sand  and  the  calculated  value. 

Acoustic  measurements  made  in  the  field  underestimated  concentration  by  20%  when 
compared  with  measurements  made  by  OBS  instruments.  This  difference  was  less  in  low 
energy  conditions.  The  near  bottom  acoustic  estimates  of  grain  size  in  the  field  were  within 
6-18%  of  the  mean  bottom  sediment  size,  and  showed  a  vertical  decrease  in  sand  size  of 
25%  over  50  cm,  which  is  consistent  with  other  field  studies. 
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III.  CALIBRATION 


A.  CASP  DESCRIPTION 

The  Coherent  Acoustic  Sediment  Probe  (CASP)  is  a  sediment  flux  sensor  which  was 
designed  and  built  by  Timothy  P.  Stanton  of  the  Naval  Postgraduate  School,  Monterey,  CA, 
in  1992/1993.  The  instrument  is  both  a  coherent  acoustic  doppler  velocity  meter  and  a 
sediment  profiler.  The  velocity  measurements  are  made  with  three  5.2  MHz  transducers 
mounted  so  that  the  beams  intersect  25  cm  in  front  of  the  instrument  (Figure  3.1).  Operation 
in  both  monostatic  and  bistatic  modes  allows  measurement  of  all  three  components  of 
velocity  at  8  mm  resolution  (Stanton  1993).  A  1.4  MHz  transducer  located  in  the  center  of 
the  instrument  and  the  three  slanted  beam  5.2  MHz  transducers  provide  a  sediment  profiling 
capability  in  1 .68  cm  bins  to  a  range  of  1  m.  Sediment  size  discrimination  can  be  inferred 
by  comparison  of  the  5.2  and  1.4  MHz  signals  at  the  intersection  point. 

Several  other  instruments  are  located  within  the  CASP  transducer  head  to  ensure 
correct  interpretation  of  the  data.  A  0.1  °  accuracy  two-axis  tilt  sensor  enables  vertical 
transformation  of  vector  quantities  into  an  earth  referenced  coordinate  system,  and  a  three 
axis  accelerometer  detects  vibrations  or  other  contaminating  movement  of  the  sensor. 
Finally,  a  high  resolution  pressure  transducer  provides  instrument  depth  and  wave  height. 

The  data  frames  output  from  the  CASP  contain  a  complex  pattern  of  both  monostatic 
and  bistatic  sampling  modes.  In  the  bistatic  mode,  a  burst  of  32  short  (0.8  cm  length  in 
seawater)  pulses  are  transmitted  by  a  single  5.2  MHz  transducer,  while  the  other  two  5.2 
MHz  transducers  “listen”  to  the  backscattered  energy  from  the  beam  intersection.  This 
continues  in  round  robin  sequence  for  each  of  the  three  5.2  MHz  transducers  (See  Stanton 
1993,  section  3  for  complete  coverage  of  velocity  measurements).  A  monostatic  sample  for 
each  of  the  four  beams  is  interleaved  within  this  sequence.  A  triplet  of  longer  ( 1 .7  cm)  bursts 
is  transmitted  by  each  of  the  four  transducers  in  sequence,  and  the  complex  backscatter 
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Figure  3.1  Schematic  of  the  CASP  acoustic  beam  geometry.  (Stanton.  1993) 
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amplitudes  are  sampled  every  1.68  cm  out  to  a  range  of  1.2  m.  The  full  sequence  of  3 
bistatic  and  3  monostatic  transmissions  is  repeated  at  a  rate  of  36  Hz.  In  the  monostatic  data 
record,  every  fourth  1.68  cm  “bin”  is  used  to  store  other  information,  so  that  the  actual 
sampling  profile  appears  as  in  Figure  3.2. 
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Figure  3.2  CASP  Bin  Map 


The  backscattered  signals  pass  through  a  high  gain,  wide  dynamic  range  amplifier  and  mixer 
which  shift  the  ultrasonic  signals  to  base  band  frequencies.  The  complex  amplitude  signals 
are  then  sampled  at  pre-programmed  delays  to  form  the  raw,  digitized,  coherently  sampled 
acoustic  data  stream.  The  resulting  processed  power  data  is  stored  as  1000  times  the  base 
ten  logarithm  of  the  counts  (voltage)  squared.  In  the  following  sections,  C  (counts)  is  used 
as  the  analogue  for  sound  pressure  ampUtude,  in  contrast  to  V  (volts)  as  used  in  the  previous 
“background”  section. 


B.  CALIBRATION  THEORY 


As  discussed  in  the  background,  inversion  of  backscatter  signals  to  sediment 
concentration  measurements  can  be  approached  either  empirically  or  through  more  analytical 
acoustic  modeling.  Because  the  sand  at  Duck,  N.C.  is  historically  homogeneous  at  the  bar, 
and  samples  are  easily  taken,  and  an  empirical  technique  is  used  for  this  calibration.  For  low 
concentrations  (approximately  10  kg'm’^),  the  dual  frequency  profiler  will  allow  some 
discrimination  of  changes  in  mean  grain  size. 
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In  the  laboratoiy,  measurements  of  known  concentrations  of  sample  sand  (from  0. 13 
to  98.9  kg  m'^)  are  made  in  the  farfield  of  the  acoustic  transducer.  The  sand  is  contained  in 
a  cylindrical  test  vessel  and  the  concentration  is  kept  homogenous  by  vigorous  mixing  (see 
laboratory  description).  Calibration  and  attenuation  coefficients  are  determined  for  each 
frequency  from  the  backscatter  levels  at  each  range. 

The  governing  equation  for  acoustic  behavior  through  a  water/sand  medium  can  be 
expressed  as: 


— F[«(a)/„(x)]— exp[-2  (3.1) 


This  equation  differs  from  Equation  2.4  in  the  inclusion  of  the  range  spreading  term  and  the 
exclusion  of  the  sinh  term.  The  within-bin  attenuation  represented  by  the  sinh  term  is  instead 
compensated  for  during  the  integration  algorithm  by  integrating  through  one  half  of  the  bin 
undergoing  calculation.  The  acoustic  response  function  F[ ]  has  a  frequency  dependent 
behavior  due  to  the  size  distribution  and  backscatter  form  function,  { x )  ,  which  was 
determined  empirically  by  Hay  (1991).  For  a  lognormal  sediment  size  distribution,  the  final 
response  function  can  be  expressed  (Hay  and  Sheng,  1992)  in  terms  of  the  geometric  mean 
radius  of  the  particles,  ,  the  acoustic  wave  number  and  the  grain  size  distribution  : 


F(X,a^)=k^ 


1/2 


-(Ina-lna 

1/-(^)|  exp[ - ^]d  \na 

i _ 

Z  -(Ina-lna  )^ 

la  exp[ - Ina 

i  21n^O 


(3.2) 


where  a=equivalent  radius  of  the  particle  and  X=kfi^  and  \fjx)  \  is  as  described  in  Equation 
2  7 
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The  response  functions  for  both  CASP  operating  frequencies  are  illustrated  for 
varying  mean  sediment  sizes  and  lognormal  distribution  widths  (a)  in  Figure  3.3.  The  sand 
sampled  from  the  bar  and  trough  area  of  Duck,  NC,  and  used  for  the  calibration  has  a  mean 
geometric  radius  of  .09  mm  and  a  standard  deviation  of  1 .2.  Note  that  for  sand  of  these 
characteristics,  the  5.3  MHz  response  is  largely  constant  with  radius,  but  can  vary  as  much 
as  10%  with  the  distribution  width.  The  1.3  MHz  beam  response  is  more  susceptible  to 
variations  in  mean  particle  size,  and  can  vary  up  to  20%  with  standard  deviation.  The 
calibration  is  made  assuming  a  constant  grain  size  distribution  with  fixed  mean  and  width 
so  that  this  form  factor  is  empirically  represented  by  the  single  frequency  calibration 
coefficient,  C„  .  Similarly,  since  the  instrument  characteristics  are  the  same  in  the 
laboratory  and  the  field,  the  instrument  sensitivity,  S,  which  depends  on  frequency, 
directivity  and  beam  characteristics,  and  on  the  form  function  of  the  object  illuminated,  will 
also  be  represented  in  the  calibration  coefficient. 

Once  these  simplifications  are  applied,  the  resulting  calibration  equation  is: 


M=C^* 


where 


(3.3  a,b) 


Af  is  the  mass  concentration  (kgm'^,  C  is  measured  counts,  C„  is  the  calibration  coefficient, 
and  A  is  the  integrated  compensation  factor  integrated  from  the  transducer  face  to  the 
measurement  bin  at  range  r. 
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Figure  3.3  Response  functions  for  1.3  MHz  and  6.3  MHz,  assuming  a  lognormal  particle 
size  distribution,  as  a  function  of  the  particle  mean  geometric  radius  for  values  of  o  =  1  1 
1.2,  1.3,  1.4,  1.5. 


C.  LABORATORY  DESCRIPTION 


1.  Tank 

The  CASP  was  immersed  in  a  600  gallon  acrylic  tank  filled  with  fresh,  filtered  water 
for  the  calibration  procedures.  A  1 6  liter  cylindrical  acrylic  test  vessel,  positioned  inside  the 
large  tank,  was  used  to  contain  controlled  masses  of  sample  sediment  (Figure  3.4).  The  test 
vessel  was  equipped  with  an  acoustically  transparent  plastic  film  window,  a  controllable 
speed  propeller  stirrer,  and  two  downward  facing  OBS’s  inserted  in  an  O  ring  sealed  sidewall 
opening.  The  vessel  was  designed  so  that  the  measured  masses  of  sand  were  homogeneously 
mixed  throughout  the  volume  with  no  level  surfaces  or  stagnation  points.  The  CASP  was 
aimed  at  the  transparent  window  so  that  the  beam  intersection  point  was  approximately  2  cm 
inside  the  test  vessel. 

2.  Sediment 

Sand  was  obtained  from  the  trough  and  bar  (top)  region  of  the  nearshore  zone  at  the 
CERC-FRF  in  Duck,  N.C.  in  both  August,  1994,  and  June,  1995.  During  the  DUCK94 
experiment,  sand  was  sampled  20  October  1 0  meters  south  of  where  the  CASP  operated 
(Stauble,  1995).  Results  of  the  grain  size  lognormal  distributions  are  summarized  in  Table 
3.1. 


Bar/Trough 

B-Oct  94 

B-Aug  94 

B-Jun  95 

T-Oct  94 

T-  Jun  95 

Mean 

.  1 42  mm 

.184  mm 

.191  mm 

.144  mm 

.192  mm 

Sigma 

1.25 

1.23 

1.2 

1.36 

1.2 

Table  3.1  Lognormal  Mean  Diameter  and  Standard  Deviation  of  Sand  from  Duck,  NC 


1.  COHERENT  ACOUSTIC  SEDIMENT  PROBE  (CASP) 

2.  ACOUSTKALLY  TRANSPARENT  MEMBRANE 

3.  ACRYLfC  NOZZLE  AND  DUCTED  PROPELLER 

4.  ACRYLIC  TEST  VESSEL 

5.  VARIABLE  SPEED  MOTOR 

6.  OPTTCAL  BACKSCATTER  SENSOR 

7-  ACRYLIC  LABORATORY  ACOUSTIC  TEST  TANK 
0.  SUSPENDED  SEDIMENT 


Figure  3.4  Schematic  of  the  test  tank,  CASP,  and  the  mixing  vessel  with  propeller  stirrer 
(Stanton.  1993). 
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The  lognormal  fit  of  the  sand  used  for  the  calibrations  is  shown  in  Figure  3.5.  The 
calibration  sand  was  sifted  to  eliminate  particles  larger  than  .85  mm  to  prevent  interference 
with  the  propeller  mechanism.  After  drying,  the  sand  was  weighed  in  roughly  sub  octave 
increments  so  that  the  concentration  in  the  test  vessel  would  increase  slowly  to  25  kg  m'\ 
then  more  rapidly  to  100  kg  m'l 

3.  Procedure 

Once  the  GASP  and  test  vessel  were  positioned,  a  dry  run  was  made  to  ensure  proper 
alignment  and  check  for  acoustic  reverberation  interference.  The  calibration  began  with  two 
cases  of  no  sediment,  one  still  and  one  stirred,  and  continued  with  incremental  additions  of 
sand.  Each  increment  was  measured  for  three  minutes,  and  one  to  two  minutes  were  taken 
between  increments  for  addition  of  sand  and  stirring.  Several  calibration  runs  were 
conducted  with  both  August  1994  and  June  1995  sand  samples. 

4.  Calibration  Method  Sensitivities 

The  consistency  of  the  calibration  data  was  affected  significantly  by  two  experimental 
variables;  (1)  the  entrainment  of  micro-bubbles  on  the  sand  particles  as  they  were  added  to 
the  test  vessel;  and  (2)  the  position  of  the  test  vessel  window  within  the  GASP  data  bin  map. 
These  considerations  wilt  be  discussed  here.  The  performance  of  the  stirrer  and  the  resulting 
homogeneity  of  the  mixture,  and  the  consistency  between  runs  using  the  two  sand  samples 
is  also  addressed. 

Rapid  addition  of  sand  to  the  test  vessel  resulted  in  an  unacceptable  degree  of  scatter 
due  to  microbubble  entrainment,  even  after  a  one  to  two  minute  settling/stirring  period 
between  recordings.  Even  after  rinsing  the  sand  with  a  chemical  wetting  agent,  the  most 
consistent  runs  were  made  when  the  sand  was  carefully  and  slowly  tapped  into  the  cylinder. 
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An  extreme  example  of  the  effect  of  this  on  the  attenuation  coefficient  is  shown  in  Figure 
3.6.  This  effect  can  also  be  observed  when  the  three  minute  recording  time  is  divided  into 
one  min  segments.  As  seen  in  Figure  3.7,  the  attenuation  coefficient  (especially  at  higher 
concentrations)  decreases  with  time,  suggesting  the  gradual  dissipation  of  bubbles  during  the 
stirring  process. 

To  examine  the  effect  of  the  window  position  within  the  bin  map,  the  CASP  was  run 
through  a  computer  controlled  grid  pattern  in  front  of  the  test  vessel  while  a  concentration 
of  10  kgm'^  was  maintained  in  the  vessel.  Figure  3.8  is  a  three  dimensional  representation 
of  the  backscatter  power  (counts^2)  response  as  the  CASP  backs  away  from  the  center  of  the 
window.  The  window  is  located  between  0.25  and  0.50  meters  as  the  CASP  retreats.  The 
rear  peak  is  the  back  of  the  test  vessel,  and  the  peaks  in  front  of  the  window  represent 
reverberation  responses.  Although  these  reverberations  do  not  affect  the  voltage  level 
received  in  the  test  vessel  segment  of  the  beam,  they  must  be  corrected  for  when  the  data  is 
used  for  the  along-beam  calibration  integration.  When  the  window  is  within  one  of  the  large 
bins  of  the  map,  the  voltage  reading  can  be  divided  between  two  bins  to  an  unknown  degree 
(Figure  3.9).  Although  this  percentage  can  be  estimated,  the  situation  causes  an  uncertainty 
factor  in  the  calculation  of  the  Cn,  coefficient.  This  will  be  discussed  further  in  the 
calculation  section. 

The  performance  of  the  stirrer  and  the  resulting  homogeneity  of  the  mixture  were  also 
investigated  as  possible  sources  of  error.  The  absolute  concentration  is  verified  by  readings 
from  two  OBS's,  mounted  one  above  the  other  at  the  level  of  the  window  and  just  below  the 
surface  of  the  water.  OBS  measurements  of  sediment  concentration  for  each  calibration  run 
agree  well  with  controlled  sediment  concentrations  (Figure  3.10).  The  one  minute  averaged 
OBS  measurements  were  consistent  within  1  kg  m‘^  below  50  kg  m'\  and  varied  less  than 
5  kgm'^  from  50  to  100  kg  m'^ ,  indicating  that  the  mixture  was  indeed  homogeneous.  If 
the  stirrer  was  offset  within  the  test  vessel,  an  inconsistent  gradient  through  the  beam  would 
result.  This  condition  could  be  detected  in  the  data  by  varying  the  number  of  bins  through 
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Figure  3.7  Evidence  of  microbubble  dissipation  with  time  for  June  calibration.  Beam  1 
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Figure  3.9  Backscatter  profile  for  1.3  MHz  during  calibration  grid  test.  A  single  peak 
indicates  that  the  first  bin  is  completely  full  of  sediment;  a  sloped  peak  occurs  when  the 
window  intersects  the  first  bin.  The  linestyles  represent  different  CASP-to-window  ranges. 


actual  concentration  (g/l) 


Figure  3.10  OBS  readings  for  June  and  August  calibration  runs. 


which  the  attenuation  coefficient  was  calculated  and  noting  inconsistencies,  but  careful 
placement  of  the  stirrer  avoided  this  situation.  The  experimental  design  and  resulting 
homogeneity  of  the  sand  volume  was  therefore  considered  sufficient  for  the  calibration 
technique. 

Any  variation  between  calibration  runs  using  each  of  the  two  sand  samples  could  be 
accounted  for  by  the  above  considerations;  minor  differences  in  the  distribution  of  the  grain 
sizes  were  not  resolved  in  the  calibration.  Two  calibration  runs,  one  using  August  sample 
bar  sand  and  the  other  June  sample  trough  sand  (hereon  referred  to  as  Aug  and  Jun  cal  runs), 
are  presented  to  illustrate  the  procedure  and  analyze  its  effectiveness. 

D  COEFFICIENT  DETERMINATION 

The  piece-wise  bin  integration  form  of  Equation  3.3b  is: 

r 

A= — exp(-2  / ci^+aslope*mass(bin-\)*dr  (3.4) 

^  D 

where  dr=l/2*length(bm-l)  +  V2*length(bm)  and  the  aslope  coefficient  is  applied  to  the 
mass  calculated  for  the  previous  bin  to  determine  the  appropriate  value  for  the 
compensation.  Two  implementations  of  the  integration  routine  are  used.  First,  an 
attenuation-only  correction  can  be  applied  to  the  counts^2  values.  In  this  case,  known 
concentration  is  used  for  the  “mass(bm-l)”  of  the  compensation  A  of  Equation  (3.4).  No  C„ 
coefficient  is  applied,  so  that  M  equates  to  corrected  counts'^2  vice  corrected  mass.  The 
second  implementation  uses  the  mass  calculated  from  the  previous  bin  as  would  be  used  in 
actual  measurements  without  known  sediment  concentrations. 

An  additional  procedure  within  the  integration  was  developed  to  correct  for  the 
situation  where  a  given  bin  spanned  the  test  vessel  acoustic  window.  This  window  bin 
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correction  was  implemented  with  a  conditional  loop  in  the  integration  routine  which 
identified  the  affected  bin  (due  to  the  Vi  bin  to  V2  bin  procedure,  the  bin  beyond  the  window 
is  affected  for  window  percentages  <.5)  and  appropriately  adjusted  the  dr  and  mass.  For  the 
corrected  mass  integration  technique,  known  mass  is  used  for  the  window  bin  compensation. 

The  data  were  also  corrected  for  system  noise.  The  backscatter  field  in  counts  for 
the  June  calibration  of  the  1 .3  MHz  beam  is  depicted  in  Figure  3.11.  The  square  root  of  the 
actual  power  was  used  to  best  identify  areas  of  reverberation  interference.  A  noise  floor  was 
calculated  by  averaging  the  values  in  the  bins  not  affected  by  reverberation  peaks,  and  the 
square  of  this  value  was  subtracted  from  the  raw  backscattered  power  record.  The  noise 
floor  averaged  4000  counts'^2.  Readings  between  the  CASP  and  the  test  vessel  face  were 
zeroed. 


1.  Aslope  Coefficient 

(L  Calculation 

Figure  3.12a  portrays  the  reading  in  counts  (volts)  along  the  CASP  1.3  MHz 
beam  while  the  vessel  concentration  was  held  at  6  kgm'^.  The  six  measurements  in  the 
center  of  the  test  volume  were  used  to  calculate  the  attenuation  coefficient,  a^.  The 
measured  change  in  backscattered  power  vs.  range  in  a  constant  concentration  follows  a 
logarithmic  decay  due  to  attenuation  of  pure  water  and  scattering  by  the  sand  particles, 
superimposed  on  a  1/r  radial  beam  spreading: 

C=C.— e where  a=a  +a  (3.5  a,b) 

0  ^  t  w  i 

Here  C  (counts)  «  acoustic  pressure,  and  total  attenuation  a„  is  the  contribution  due  to  fresh 
water  absorption  and  particle  scattering,  respectively. 
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Therefore, 


a  = 


.,Cr 

InC- - ) 

-lt\r 


-a 


w 


(3.6) 


and  a^is  calculated  for  each  concentration  as  the  rate  of  the  decay  of  the  range-corrected 
natural  logarithm  of  counts  over  range  (Figure  3. 12b).  The  value  of  each  is  then  plotted 
versus  concentration,  and  the  calibration  attenuation  coefficient  as/qpe(/(m*kgm'^))  is 
determined  by  linear  regression.  To  integrate  a  sediment  measurement  profile  using 
Equation  (3.4) ,  the  aslope  coefficient  is  applied  to  the  concentration  of  the  previous  bin  to 
determine  the  a^.  value,  and  the  a„  and  1/r  spreading  corrections  are  reapplied  to  calculate 
the  compensation  factor  for  that  bin. 

The  a„  value  does  not  vary  significantly  between  seawater  and  freshwater  for 
either  frequency,  but  does  change  with  frequency  and  temperature.  Values  were  calculated 
from  Fisher  and  Simmons  (1988)  and  are  shown  in  Table  3.2: 


Temp/Freq 

1.3  MHz 

5.3  MHz 

15  C 

0.0503 

0.8369 

20C 

0.0414 

0.6884 

23C 

0.0371 

0.6166 

Table  3-2.  cq.  Values  for  Frequency  and  Temperature  Variations 

b.  1.3  MHz  Beam  Attenuation 


Attenuation  coefficients  for  the  1 .3  MHz  beam  for  the  full  range  of  concentrations 
for  the  June  and  August  calibrations  are  shown  in  Figure  3.13.  A  common  characteristic 
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of  all  calibration  runs  was  the  nonlinearity  of  values  above  33  kg  m'^.  Further  analysis  of 
corrected  sediment  profiles  confirmed  a  breakdown  of  predictable  behavior  at  this  point, 
indicating  the  onset  of  multiple  scattering.  The  first  three  anomalous  points  are  due  to 
residual  background  particles  at  the  0  and  0.13  kg  m'^  concentrations.  The  ealibration 
calculations  therefore  used  concentrations  from  .5  to  25  kgm’^  (Figure  3.14).  The  aslope 
values  differed  slightly  (<0.005  /(mkgm’^))  when  four  to  six  bins  beyond  the  window  were 
used  for  the  calculation;  therefore  the  mean  of  these  values  is  used  for  the  calibration  (Table 
3.3).  Other  calibration  runs  with  minimal  micro-bubble  contamination  yielded  values  within 
.005  /(m  kg  m'^)  of  this  value. 

The  attenuation-only  integration  technique  was  applied  to  the  counts^2  calibration 
data  to  evaluate  the  performance  of  the  calculated  aslope  coefficient.  Figures  (3.15a-g) 
illustrate  this  correction  for  each  concentration.  At  bin  18(13  cm  within  the  volume),  the 
attenuation  correction  equates  to  a  60  to  80%  correction  in  mass  from  the  lightest  to  heaviest 
(25  kg m’^)  concentrations.  The  integration  performs  well  through  the  entire  1 3- 1 5  cm  range; 
standard  deviation  along  the  line  averaged  less  than  0.1  (counts^2/le8)  up  to  15  kg  m  ^  and 
less  than  0.3  (counts^2/le8)  up  to  25  kg  m‘^.  The  instability  of  the  correction  at  33  kg  m'^ 
attests  to  the  onset  of  multiple  scattering  effects. 


Cal  /  Beam 

1 

2 

3 

4 

June 

-0.0588 

-1.113 

-1.147 

-1.117 

August 

-0.0574 

-1.123 

-1.197 

-1.161 

Calibration 

-0.058 

-1.14 

-1.14 

-1.14 

Table  3.3.  Aslope  Coefficients 
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JUN  cal  Beam  1 


bin  bin 


Figure  3.15a  1.3  MHz  backscatter  voltage  corrected  for  attenuation;  o^^uncorrected 

counts'^2.  *=corrected  counts'^2.  The  two  0  g/1  examples  are  for  an  empty  test  vessel  first 
unstirred  and  then  stirred. 


c.  5.3  MHz  Beam  Attenuation 


The  effectiveness  of  the  5.3  MHz  beams  suffers  from  multiple  scattering  at  a 
much  lower  particle  concentration  (Figure  3.16);  the  aslope  calculations  were  therefore 
limited  to  8  kg•m■^  The  aslope  values  in  Table  3.3  are  again  the  mean  of  the  values 
calculated  using  four  to  six  bins  beyond  the  window  (these  values  varied  an  average  of 
.03/(mkg-m'^)  from  one  another).  The  attenuation  integration  was  conducted  with  the  partial 
window  modification;  Figures  3. 17a-d  show  the  results  of  this  correction  through  the  beam. 
The  5.3  MHz  attenuation  correction  is  consistent  only  3  bins  (5-7  cm)  into  the  medium  once 
concentrations  exceed  4  kg  m'^. 

2.  Cm  Coefficient 

a.  Calculation 

The  calibration  coefficient  (Equation  3.3a),  quantifies  the  linear  proportionality 
of  mass  to  the  square  of  pressure  amplitude  counts  measured  by  the  instrument  (corrected 
by  spreading  and  attenuation.).  Several  methods  exist  to  calculate  this  proportionality  and 
include  direct  fit,  attenuation-corrected  fit,  and  optimization  fit  techniques. 

First,  and  most  simply,  is  the  direct  linear  fit  of  the  counts^!  reading  in  the  first 
bin  of  the  test  vessel  to  the  known  concentration.  As  this  fit  is  made  to  increasing  maximum 
concentrations  (4  to  25  kgm'^),  its  value  changes  slightly,  a  behavior  attributed  to  the 
approach  of  the  multiple  scattering  regime.  Ideally,  the  window  of  the  test  vessel  is  located 
at  the  very  beginning  of  a  data  bin,  so  that  the  first  bin  for  calculation  is  represented  by  a  full 
bin  of  sediment.  In  this  case,  the  values  for  this  first  bin  need  only  be  corrected  for  the  water 
attenuation  to  the  window,  and  then  fit  to  the  known  concentration.  The  experimental 
technique  emphasized  correct  positioning  of  the  1.3  MHz  beam  window  bin,  which 
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Cts'^2(/1e8)  cts''2(/1e8) 


cts''2(/1e8)  cts^2(/1e8) 


JUN  cal  Beam  2 


5.1g/l  6,1g/l 


Figure  3.17c  5.3  MHz  backscatter  voltage  for  Beam  2  corrected  for  attenuation; 

o=uncorrected  counts^2,  *=corrected  counts'^2. 


facilitated  the  1.3  MHz  calculations,  but  because  of  the  geometry  of  the  instrument,  resulted 
in  a  50  to  80%  window  position  for  each  of  the  5,3  MHz  beams  (Figure  3.18a).  An 
estimate  of  this  window  percentage  was  made  by  extrapolating  the  a,  log-linear  attenuation 
fit  back  to  the  window  bin  for  each  concentration  (Figure  3 . 1 8b),  and  comparing  this  value 
(in  counts  space)  with  the  actual  counts  reading  in  the  window  bin  (Figure  3.18c).  This 
extrapolated  window  bin  value  was  also  examined  for  use  in  the  direct  calculation  method. 
Some  error  exists  because  the  attenuation  experienced  in  the  first  partial  bin  lowers  the 
counts  in  the  subsequent  bins  from  which  the  extrapolation  is  performed;  therefore  the 
correction  is  not  a  completely  true  representation  of  absolute  counts  read  by  the  instrument 
at  higher  concentrations. 

The  attenuation-corrected  fit  calibrates  known  mass  with  the  attenuation-corrected 
counts^2  values  described  above.  The  calibration  coefficient,  again  a  linear  regression  to 
concentration,  can  be  calculated  using  either  a  single  value  of  the  window  bin  or  one  bin 
beyond  the  window,  or  an  average  of  the  values  in  a  selected  number  of  bins  beyond  the 
window.  The  main  advantage  to  this  technique  is  its  ability  to  correct  for  within-bin  and 
partial  window  attenuation. 

The  optimization  fit  technique  applies  a  range  of  C„  coefficients  to  the  mass 
calculation  in  order  to  determine  which  value  best  corrects  the  data  to  the  known  sediment 
concentrations.  The  problem  is  two-dimensional;  the  performance  of  the  coefficient  must 
be  satisfactory  for  the  range  of  concentration  at  the  first  bin,  and  at  bins  further  within  the 
test  volume.  Therefore  the  optimization  criteria  selected  is  the  average  least  squares  fit  of 
corrected  to  actual  concentration  for  the  three  bins  beyond  the  window.  Results  are 
described  below. 
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b.  1.3  MHz  Mass 


Table  3-4  lists  the  C„  values  calculated  from  the  three  methods.  The  window-bin 
counts'^2  values  used  for  the  direct  fit  calculation  are  shown  in  Figure  3,19.  The  window 
bin  for  the  1.3  MHz  calibration  runs  was  examined  for  partial  window  effects  with  the 
back-extrapolation  method  described  previously.  At  concentrations  above  5  kgm'^, 
calculations  indicate  that  the  instrument  was  seeing  a  full  window  of  sediment.  Below  5 
kgm'^,  however,  the  reading  is  greater  than  the  expected  value.  This  effect  is  thought  to  be 
due  to  side  lobe  or  small  scale  reverberation  interference,  and  would  cause  the 
calculated  C„  coefficient  to  be  less  than  its  actual  value,  especially  when  the  fit  is  done  only 
to  these  low  concentrations. 

The  attenuation-corrected  fit  method  used  the  counts^2  value  one  bin  beyond  the 
window,  and  the  average  of  the  values  in  four  bins  beyond  the  window  (Figure  3.20). 
Because  of  the  consistency  of  the  attenuation  con  ection  through  the  medium,  the  averaging 
technique  with  more  data  points  is  favored. 

The  C„  values  from  the  optimization  technique  were  in  close  agreement  and 
values  ranged  from  3.6  to  3.9  for  fits  from  4  to  25  kgm'^.  The  corrected  mass  accuracy  for 
the  three  bins  beyond  the  window  (approximately  7  cm)  averaged  8%  for  the  calculation  to 
8  kg^m'^  and  1 1%  for  the  calculation  to  25  kgm'^  (Figure  3.21).  The  appropriate  correction 
was  therefore  chosen  to  be  3.7  (x  e-8  /kgm'^);  the  corrected  surfaces  for  both  calibrations 
using  this  value  are  presented  in  Figures  3.22  and  3.23.  Application  of  this  coefficient  to 
the  noise  floor  backscatter  level  produced  a  minimum  measurable  concentration  of  O 
(1x10’''  kg  m'^). 
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concentration  (g/l) 


Figure  3.20  Attenuation-Correction  method  backscatter  values  for  1.3  MHz 


corrected  mass  (g/l) 


corrected  mass  (g/l) 


Figure  3.23  Corrected  mass  for  Beam  1  August  run  using  calibration  C,„. 
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Method: 

Direct 

Att-Corr  (bin) 

Att-Corr  (avg) 

Optimization 

Limit: 

8g/I 

25g/I 

8g/l 

25g/l 

8g/I 

25  g/i 

8g/l 

25g/I 

June 

2.7 

3.1 

4.1 

4.4 

3.3 

3.8 

3.6 

3.9 

August 

2.8 

3.2 

3.8 

4.2 

3.3 

3.6 

3.5 

3.9 

Table  3.4  Beam  1  C„  Coefficient  Calculated  Values 
c.  5.3  MHz  Mass 

C„  values  for  the  5.3  MHz  calculations  are  given  in  Table  3.5.  The  procedures 
for  the  5.3  MHz  beams  were  complicated  by  the  partial  window  geometry  of  the  experiment 
configuration,  as  described  above.  Figure  3.24  displays  the  window  bin,  bin  beyond  the 
window,  and  corrected  window  bin  readings  for  beam  2  of  the  June  calibration.  Direct  fits 
done  to  the  corrected  data  yielded  values  less  than  those  of  the  optimization  method,  which 
would  follow  from  the  inaccuracy  of  the  extrapolation  method  which  was  discussed  above. 

Attenuation-correction  fit  method  calculations  were  sensitive  to  the  partial 
window  value,  and  to  the  choice  of  single  or  averaged  bin  technique.  The  results  of  this 
computation  are  illustrated  in  Figure  3 .25.  The  3-bin  average  method  yielded  a  coefficient 
1/10  lower  than  the  single  bin  beyond  the  window  calculation.  Examination  of  the  counts^2 
corrections  (Figures  3.17a-d  )  suggests  that,  due  to  the  spread  of  the  corrected  values,  the 
single  bin  method  is  more  reliable.  Values  in  Table  3.5  were  calculated  using  the  single  bin 
method  and  the  window  percentage  calculated  from  the  back-extrapolation  method.  An 
inaccuracy  of  this  technique  which  may  have  more  effect  on  the  5.3  MHz  than  the  1.3  MHz 
instrument  is  the  introduction  of  the  known  mass  at  the  window  bin;  if  this  value  is  higher 
than  the  concentration  actually  seen  at  the  window  interface,  the  difference  could  cause  a 
low  C„  coefficient. 


58 


The  optimization  method  proved  to  be  the  most  reliable  technique  for  deriving  an 
accurate  C„  for  the  5.3  MHz  beams.  The  method  was  used  with  a  window  correction,  and 
best  fits  made  to  4  kgm'^  (Table  3.5).  Best  fits  to  8  kgm'^  varied  only  .01  (x  e-8  /kgm*^) 
from  these  values.  The  corrections  were  consistent  to  3  bins  beyond  the  window  (Figure 
3.26).  Both  the  June  and  August  calibrations  were  accurate  to  7-8%  of  the  actual 
concentration  to  8  kgm'^  The  values  selected  as  the  calibration  coefficients  were  therefore 
an  average  of  these  coefficients.  The  coefficients  for  the  different  beams  differ  due  to 
transducer  and  electronic-specific  gain  characteristics;  beam  4  was  consistently  lower  than 
the  other  two  beams  in  all  calibration  derivations.  A  representative  corrected  surface  for 
the  5.3  MHz  calibration  is  presented  in  Figure  3.27.  The  minimum  measurable 
concentration  for  the  5.3  MHz  beams  is  0(1x10'^  kgm'^). 


Method: 

Direct 

Att-Corr 

Optimum 

Calibration 

%Wbin 

B2  June 

0.80 

0.78 

0.95 

0.93 

0.62 

B2  August 

0.96 

0.84 

0.90 

0.93 

0.82 

B3  June 

0.86 

0.81 

1.00 

0.98 

0.58 

B3  August 

0.87 

0.88 

0.97 

0.98 

0.73 

B4  June 

0.67 

0.64 

0.78 

0.77 

0.66 

B4  August 

0.69 

0.68 

0.76 

0.77 

0.74 

Table  3.5  5.3  MHz  Beam  C„  Coefficient  Calculated  Values  to  4  kgm'^ 

3.  Full  Optimization 

The  successful  use  of  the  optimization  routine  leads  to  the  speculation  that  a 
double-optimization,  which  also  includes  likely  values  of  aslope,  would  confirm  the 
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corrected  mass  (g/l) 


range  (m) 


Figure  3.27  Corrected  mass  for  Beam  3,  August  run  using  calibration  C,„. 
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calibration.  For  the  1 .3  MHz  beam,  this  appears  to  be  so,  (Figure  3.28),  although  the  spread 
of  satisfactory  aslope  coefficients  is  quite  broad.  On  the  other  hand,  the  optimization  map 
for  the  5.3  MHz  beam  is  a  multi-valued  surface  (Figure  3.29).  This  is  attributed  to  the 
nonlinear  form  of  the  integration  (Equation  3.3),  and  the  high  attenuation  values  for  the 
higher  frequency  beam.  The  optimization  evaluates  the  fit  of  a  data  line  which  originally 
curves  downward  toward  higher  concentration  values;  the  coefficient  determines  the 
slope  of  the  corrected  line,  and  the  aslope  coefficient  determines  the  degree  to  which  the  data 
points  are  brought  up  to  that  line.  It  is  therefore  conceivable  that  more  than  one  set  of  values 
would  optimize  the  correction,  especially  in  the  heavily  attenuated  5.3  MHz  case.  Because 
the  attenuation  correction  described  here  is  considered  appropriate,  and  the  precise 
optimization  values  depend  upon  the  criteria  chosen  for  optimization  (averaged  number  of 
bins  beyond  the  window,  etc),  the  double  optimization  is  not  used  for  attenuation  coefficient 
calculation. 

E.  GRAIN  SIZE  DISCRIMINATION 

The  grain  size  discrimination  capability  of  the  CASP  is  based  on  the  implicit 
inclusion  of  the  form  factor  F[  ]  in  the  calibration  coefficient  C„.  Because  both  sonar 
frequencies  were  calibrated  to  the  same  sand  sample,  the  ratio  of  the  corrected  masses  should 
be  unity.  If  we  assume  that  the  sand  size  spectrum  is  represented  by  a  lognormal 
distribution,  the  absolute  ratio  of  the  actual  form  factors,  (Figure  3.30),  would  then  be 
normalized  by  a  secondary  calibration  coefficient,  (in  this  case  equal  to  0.625)  so  that 
the  ratio  at  0.09  mm  radius,  1.2  standard  deviation  (i.e.  calibration  sand)  would  be  one. 
Should  the  ratio  in  the  field  deviate  from  unity,  this  secondary  coefficient  would  be  reapplied 
to  estimate  the  actual  grain  size.  It  follows  from  the  above  calibration  description  that  the 
performance  of  the  grain  size  discriminator  will  be  limited  by  the  5.3  MHz  beam.  Figures 
3.31  and  3.32  show  that  this  is  the  case.  For  both  calibration  runs,  the  mass  ratio  is 


64 


3 


0.5 


0.4 


0.6  0.8  1  1.2  1.4  1.6 

Mean  Geometric  Radius  (m) 


Figure  3.30  Theoretical  values  of  F(X5,Og)/F(X|,ag)  as  a  function  of  ag  for  Og  =  1 . 1 ,  1 .2, 
1.3,  1.4,  1.5. 
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F5/F1 


consistent  to  approximately  10%  at  concentrations  less  than  8  kg  m'^  and  distances  less  than 
7  cm.  It  is  reasonable  to  expecct  that  the  instrument  will  be  able  to  detect  mean  grain  size 
varioations  of  more  than  approximately  0.03  mm  in  light  sediment  conditions. 
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IV.  FIELD  DEPLOYMENT 


A.  SITE  DESCRIPTION 

The  CASP  was  deployed  as  part  of  a  multi-instrument  package  installed  on  a 
movable  sled  during  the  DUCK  94  experiment  in  October  1994  at  the  U.S.  Army  Corps  of 
Engineers’  Field  Research  Facility  (FRF)  in  Duck  NC.  Dr.  E.B.  Thornton  and  T.P.  Stanton 
were  co-principal  investigators  for  the  instrument  sled.  The  CASP  was  mounted  on  a 
computer  controlled  hydraulic  arm  which  allowed  it  to  be  positioned  at  selected  heights  from 
the  bottom.  The  sled  was  deployed  daily  from  6  to  20  Oct  along  a  cross-shore  transect 
starting  at  or  seaward  of  the  sand  bar  and  towed  across  the  surf  zone.  Data  were  acquired  by 
stopping  at  4  to  8  stations  along  the  transect  each  day  as  the  sled  was  dragged  shoreward. 
At  each  station,  the  CASP  performed  an  approximately  one  hour  profile  which  began  by 
positioning  the  beam  intersection  point  (25  cm  from  the  face  of  the  CASP)  as  closely  as 
possible  to  the  bottom,  and  stepping  the  instrument  up  in  approximately  one  bin  (1 .68  cm) 
increments.  Each  increment,  termed  a  substation,  was  two  to  ten  minutes  long.  At  each 
station,  five  to  twenty  substations  were  recorded. 

Data  records  from  Oct  20th  were  chosen  to  illustrate  CASP  capabilities  because  of 
its  strong  low  frequency  wave  forcing  and  absence  of  bubble  contamination.  Of  four  stations 
recorded  that  day,  the  farthest  station  from  the  beach.  Station  1,  contained  the  most 
consistent  substation  records.  The  positions  of  the  four  stations  are  noted  on  a  bathymetric 
profile  of  the  experiment  area,  Figure  4.1.  At  Station  1  the  sled  was  located  on  the 
shoreward  slope  of  the  offshore  bar,  on  an  incline  of  approximately  2  cm/m.  The  small  scale 
bathymetry  of  the  area,  measured  with  altimeter  and  side  scan  sonar  mounted  on  a 
specialized  nearshore  survey  vehicle  (CRAB),  was  relatively  smooth,  and  consisted  of  one 
to  two  cm  ripples  approximately  every  five  meters.  Grain  size  analysis  performed  1 0  meters 
south  of  the  sled  transect  (Stauble,  1995)  revealed  a  sand  population  with  a  slightly  smaller 
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geometric  mean  radius  than  the  calibration  sand  (Table  3.1).  The  narrow  band  wave  field 
on  October  20th  had  an  incidence  angle  of  4“,  with  a  regular  period  of  14  seconds,  group 
period  of  approximately  2  minutes,  and  significant  wave  height  of  one  meter  at  8  meters 
depth.  The  wave  frequency-direction  spectrum  for  that  day.  Figure  4.2,  illustrates  the 
narrow  banded  directional  spectrum.  The  longshore  current  was  negligible  and  varied  by 
±  0.2  m/s. 

B.  SEDIMENT  VS.  VELOCITY  PHASE  ANALYSIS 

1.  Background  Theory 

An  analysis  of  the  correlation  of  sediment  and  cross-shore  flow  with  height  above 
the  bed  is  presented  for  day  20  to  illustrate  the  capabilities  of  the  CASP.  Current  theory  for 
suspended  sediment  distribution  in  pure  oscillatory  flow  over  a  smooth  bed  hypothesizes  that 
the  maximum  concentration  just  above  the  bed  occurs  simultaneously  with  the  maximum 
bed  shear  stress.  The  time  of  this  maximum  bed  shear  stress  has  been  calculated  as  a  phase 
lag  from  the  maximum  wave  velocity  outside  the  boundary  layer  [Fredsoe  and  Diegaard, 
1992].  Above  the  bed,  maximum  concentration  lags  the  maximum  bed  shear  stress  to  an 
increasing  degree  with  height.  Figure  4.3a  [Justesen  and  Fredsoe  1985]  shows  the  expected 
concentration  relative  to  bed  concentration  at  1.8  cm  above  bed  level  for  a  full  wave  cycle 
as  predicted  by  three  hydrodynamic  approaches  described  in  Fredsoe  and  Diegaard  [1992]. 
Each  of  these  models  assumes  forcing  by  a  symmetrical  sinusoidal  wave  form,  and  predicts 
that  the  maximum  concentration  will  be  equal  under  both  inshore  and  offshore  flow 
conditions.  Variation  of  concentration  with  wave  phase  at  different  heights  above  the  bed 
was  calculated  by  Fredsoe  et  al  [1985]  (Figure  4.3b)  using  time- variant  eddy  viscosity  and 
bed  concentration  based  on  the  Shields  parameter  [Fredsoe  and  Diegaard,  1992],  Using  a 
grain  diameter  of  0.18  mm,  the  Z/K„  (depth1?ed  roughness)  parameter  equates  to  the 
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Figure  4.3  a)  Comparison  between  three  hydrodynamic  approaches  to  estimating  eddy 
viscosity.  (From  Justesen  and  Fredsoe,  1985).  b)  Variation  in  concentration  with  time  a 
different  heights  from  the  bed.  (From  Fredsoe  et  al.,  1985). 
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following  heights  (Table  4. 1): 


Z/Kn 

0.800 

1.5 

2.7 

10.5 

23.4 

37.9 

Z  (cm) 

.036 

.067 

.121 

Al 

1.05 

1.70 

Table  4. 1  Z  for  Z/K„  using  grain  diameter  of  0. 1 8  mm 


With  the  addition  of  time  variant  eddy  viscosity,  the  maximum  concentration  at  1.7  cm 
shifts  to  approximately  7t/2,  or  time  of  maximum  cross-shore  velocity. 

Sediment  behavior  over  a  rippled  bed  becomes  much  more  complex,  and  requires 
numerical  solution  of  the  flow  equations.  Using  a  discrete  vortex  model  under  symmetrical 
forcing,  Hansen  et  al  [1991, 1994]  calculated  the  flow  field  over  a  ripple  bed  at  various  wave 
phases.  Then,  with  a  Lagrangian  formulation,  the  paths  of  the  individual  sand  grains  were 
followed.  The  calculated  pattern  of  suspended  particles  at  two  phases  of  a  wave  cycle  are 
illustrated  in  Figure  4.4. 

2.  Phase  Mapping  Method 

A  typical  time  series  of  cross-shore  velocity  and  sediment  concentrations  at  one  to 
three  bins  (2  to  7  cm)  above  the  bottom  at  Station  1  is  shown  in  Figure  4.5  (positive 
velocities  are  offshore).  The  correlation  of  maximum  onshore  velocity  (therefore  maximum 
bed  stress)  with  concentration  at  the  bottom  (bin  1),  and  the  decorrelation  with  height  (bins 
2  and  3  coirespond  to  3.4  and  5.1  cm  above  the  bed),  is  readily  apparent.  In  order  to  more 
quantitatively  analyze  this  correlation,  a  phase  mapping  display  is  used.  The  velocity  time 
series  was  broken  into  wavelengths  using  a  zero-down  crossing  technique,  so  that  a  full 
period  is  defined  as  the  beginning  of  onshore  flow  (negative  v)  to  the  end  (post-maximum) 
of  offshore  flow.  The  concentration  (or  u,v,w  magnitude)  time  series  can  then  be  separated 
into  these  same  time  defined  full  periods,  and  the  data  ensemble  averaged  relative  to  phase. 
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Figure  4.4  Modeled  path  of  suspended  sediment  particles  over  a  rippled  bed  at  two  phases 
of  a  wave  period.  (From  Hansen  et  al.,  1991). 
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A  typical  1 0  minute  substation  record  consists  of  over  40  cycles  for  the  14  second  period 
waves  observed  on  the  20th  of  October.  For  this  study,  each  360  degree  period  was  divided 
into  5  degree  segments  for  averaging.  The  correlation  of  the  cross-shore  velocity  magnitude 
with  its  own  phase  was  first  examined  to  determine  the  shape  of  the  cross-shore  wave,  and 
to  detect  any  variation  of  the  velocity  with  elevation  from  the  bed.  At  Station  1,  six 
substation  measurements  were  taken  starting  with  velocity  measurement  at  the  bed  (Bin  12, 
25  cm  from  the  face  of  the  CASP),  and  at  elevations  of  up  to  12  cm  from  the  bed.  The  shape 
of  the  velocity  wave  was  consistent  at  each  height  (Figure  4.6).  It  is  uniformly  asymmetrical, 
with  the  offshore  maximum  at  56  degrees  ±3  degrees  through  the  water  column.  The 
assymmetry  of  these  narrow  banded  waves  are  indicative  of  strong  nonlinearities  in  shallow 
water  waves.  The  velocity  magnitudes  at  the  bed  were  somewhat  less  than  those  above,  as 
expected. 


3.  Concentration  Data 

Results  of  the  mapping  of  mass  concentration  as  a  function  of  wave  phase  are  shown 
for  four  substations  in  Figures  4.7  to  4.10.  At  the  bed,  the  maximum  mass  concentration  is 
cleanly  correlated  with  the  maximum  onshore  velocity,  but  consistently  decorrelates  one  bin 
(3-4  cm)  from  the  bed.  Note  that  the  concentration  ranges  span  from  0(50)  to  0(  1 0"^)  kg  m‘ 
The  behavior  observed  in  these  plots  can  be  related  to  a  simple  conceptual  model  in  which 
sand  is  injected  into  the  water  column  from  the  boundary  layer  as  described  by  time- variant 
eddy  viscosity  theory  (Figure  4.3b),  except  that,  under  these  assymetric  waves,  the  injection 
only  occurs  at  the  onshore  phase.  Above  the  boundary  layer,  the  sand  distribution  is 
controlled  by  advection  and  turbulence  in  the  water  column. 
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Figure  4.7  Log  mass  profiles  of  phase-mapped  sediment  from  20  Oct,  Stat 
4.  The  upper  profile  is  the  bottom  bin,  with  range  from  the  bed  increasing 
Inp  mass  in  one  bin  steps  to  an  ultimate  elevation  of  1 8  cm. 
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Figure  4. 1 0  Log  mass  profiles  of  phase-mapped  sediment  fr( 
7.  The  upper  profile  is  the  bottom  bin,  with  range  from  the 
log  mass  in  one  bin  steps  to  an  ultimate  elevation  of  1 8  cm 


4.  Sediment  Flux 


Sediment  flux,  the  product  of  velocity  and  concentration,  and  can  be  calculated  with 
elevation  as  the  mean  flux  for  a  given  time  period  (here,  a  substation),  or  represented  in 
phase  space.  Mean  fluxes  for  four  substations  of  Station  1  are  presented  in  Figure  4.11. 
For  the  phase  space  calculations,  the  difference  in  the  velocity  profile  taken  at  the  bottom 
bin  and  the  profiles  taken  above  the  boundary  layer  (Figure  4.6),  can  be  accounted  for  by 
using  the  boundary  layer  velocity  profile  to  determine  sediment  flux  at  the  bottom  bin  of 
each  of  the  other  stations.  Figure  4.12  is  an  example  of  this  calculation  for  substation  5. 
The  majority  of  the  transport  occurs  in  the  bottom  boundary  layer,  with  a  small  amount  of 
onshore  transport  higher  in  the  water  column. 

C.  GRAIN  SIZE  DISCRIMINATION 

As  described  by  equation  3. 1,  the  F5/F1  value  is  proportional  to  the  square  root  of 
the  ratio  of  the  mass  measured  by  the  1 .3  MHz  sonar  to  the  mass  measured  by  the  5.3  MHz 
sonar.  The  proportionality  coefficient,  0.65,  derived  in  the  calibration  section,  was  applied 
to  the  calculated  ratio  so  that  a  corrected  F5/F1  value  of  1.6  indicates  sediment  of  the  same 
lognormal  characteristics  as  the  calibration  sediment.  Table  4.2  lists  the  corrected  F5/F1 
ratios  for  each  substation  of  Station  1,  the  height  of  the  measurement  above  the  bottom,  and 
the  max  and  mean  sediment  measured  by  the  1.3  MHz  beam.  Table  4.3  lists  these  same 
values  for  Station  2  data.  When  the  intersection  point  is  at  the  sediment  bed,  the  magnitude 
of  sediment  concentration  exceeds  the  expected  limits  for  satisfactory  performance  by  the 
grain  size  discriminator.  The  actual  depth  of  the  bottom  bin  within  the  boundary  layer  is 
also  subject  to  an  uncertainty  of  1.68  cm,  and  is  probably  responsible  for  the  difference 
between  the  ratios  at  substations  4  and  10  of  Station  1.  For  these  reasons,  the  ratios  at  the 
bottom  bin  are  not  considered  reliable. 
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Figure  4.1 1  Mean  sediment  flux  measurements  at  four  substations  of  Station  1,  Day  20 
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Substa 

Height 

(cm) 

Max 

m 

Mean 

(g/1) 

Beam  2 
F5/F1 

Beam  3 
F5/F1 

Beam  4 
F5/F1 

4 

0 

125.40 

38.0 

6.00 

9.31 

8.96 

10 

0 

141.80 

39.4 

2.95 

2.47 

2.53 

3 

3.4 

1.65 

0.16 

1,06 

1.09 

1.03 

5 

5.0 

2.52 

0.11 

1.01 

1.03 

0.97 

6 

10.0 

1.59 

0.11 

0,99 

1.00 

0.95 

7 

12.0 

1.15 

0.04 

0.98 

1.00 

0.95 

Table  4.2  F5/F1  ratios  for  Day  20,  Station  1 


Substa 

Height 

(cm) 

Max 

(g/1) 

Mean 

(g/I) 

Beam  2 

F5/F1 

Beam  3 

F5/F1 

Beam  4 

F5/F1 

2 

0 

127.10 

18.01 

927.01 

2.40 

2.38 

6 

3.4 

15.72 

1.04 

1.88 

2.03 

1.90 

3 

5.0 

1.37 

0.20 

1.96 

1.94 

1.82 

4 

10.0 

1.52 

0.38 

2.38 

2.37 

2.23 

Table  4.3  F5/F1  ratios  for  Day  20,  Station  2 

The  F5/F1  values  above  the  bed  at  Stations  1  and  2,  where  sediment  conditions  are 
well  within  those  acceptable  for  discriminator  performance,  are  notable  for  their  consistency 
per  station.  Interpreted  with  Figure  3.30,  these  values  indicate  that  the  sediment  seen  at 
Station  1  is  of  a  larger  mean  radius  (>.  1  mm),  and  the  suspended  sediment  at  Station  2  is  of 
a  smaller  mean  radius  (<.05  mm),  than  that  of  the  sand  used  for  instrument  calibration.  This 
may  be  due  to  the  higher  wave  energy  level  at  Station  1  (bar)  vs.  Station  2  (trough),  but 
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further  investigation  is  required.  Sediment  at  Station  2  also  shows  a  general  decrease  in  size 
with  height,  which  is  expected.  Bubble  entrainment  differences  would  effect  measurements 
from  station  to  station,  and  work  at  the  Naval  Postgraduate  School  continues  towards 
identifying  and  quantifying  these  effects. 
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V.  CONCLUSIONS 


An  empirical  calibration  of  a  Coherent  Acoustic  Sediment  Profiler  (CASP)  has  been 
performed  with  sediment  sampled  from  the  nearshore  zone  at  Duck,  N.C.  The  sediment  is 
characterized  by  a  lognormal  distribution  with  a  mean  grain  size  diameter  of  0.18  mm,  and 
standard  deviation  of  1.2.  The  1.3  MHz  sonar  measures  sediment  concentrations  of  0(1x10'^ 
kgm‘^)  to  25  kgm'^  at  10%  accuracy.  The  5.3  MHz  sonar  is  similarly  accurate  for 
concentrations  of  ©(IxlO'^  kgm'^)  to  8  kgm  l  The  high  concentration  limits,  25  kgm'^  and  8 
kgm agree  with  previous  semi-empirical  work  performed  using  similar  frequencies  and 
slightly  smaller  sand  size  (Sheng  and  Hay,  1988),  and  result  from  the  onset  of  multiple 
scattering  at  the  respective  frequency.  Estimation  of  the  F5/F1  form  factor  ratio  can  be  made 
by  comparison  of  concentration  measurements  by  both  frequencies,  and  provides  a  grain  size 
discrimination  capability.  This  capability  is  accurate  to  within  ±5%  of  the  F5/F1  ratio  within 
the  concentration  and  range  limits  of  the  5.3  MHz  beam. 

The  data  storage  bin  map  design  and  the  effect  of  the  partial  inclusion  of  the  test  vessel 
window  in  the  first  bin  for  the  5.3  MHz  beams  were  complicating  factors  for  calibration.  The 
bin  map  is  being  redesigned  to  minimize  the  number  of  blank  storage  bins  in  the  string.  To 
avoid  the  placement  of  the  test  vessel  window  partially  within  a  bin,  each  beam  should  be 
calibrated  separately.  A  real  time  analysis  program  has  been  implemented  in  the  laboratory 
to  allow  more  accurate  positioning  of  the  window  bin  for  any  beam  with  respect  to  the  test 
vessel  window. 

Data  from  a  day  characterized  by  narrow  banded  waves  and  negligible  longshore  current 
during  the  DUCK94  experiment  show  that  episodes  of  high  sediment  concentration  correlate 
well  with  the  maximum  onshore  velocities  within  2  cm  of  the  bottom,  but  decorrelate  above  this 
level.  A  phase  mapping  technique  is  used  to  create  ensemble  averages  of  velocity  and 
concentration  magnitudes.  Four  of  these  ensemble  sets,  consisting  of  values  from  30  to  88 
wave  periods,  are  presented  from  a  data  collection  station  located  on  a  smooth  bottom  on  the 
shoreward  slope  of  the  offshore  bar.  All  show  the  same  strong  correlation  at  the  bottom  bin. 
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and  slightly  different  sediment  patterns  higher  in  the  water  column.  Sediment  flux  calculations 
reveal  strong  shoreward  transport  at  the  boundary  layer,  and  weaker  seaward  transport  higher 
in  the  water  column.  The  grain  size  discrimination  indicates  that  the  suspended  populations 
differ  at  two  of  the  data  collection  stations. 

The  CASP  pssesses  additional  capabilites  to  allow  a  more  exhaustive  analysis  of  the 
processes  responsible  for  the  sediment  behavior  presented  here.  Vertical  velocities  measured 
by  the  CASP  are,  as  of  this  writing,  corrected  for  any  tilt  of  the  CASP  itself,  but  not  for  large 
and  small  scale  bathymetry  on  which  the  sled  might  be  resting.  Currently,  these  measurements 
provide  input  for  conceptual  modeling,  and  will  soon  be  available  for  more  quantitative  uses. 

Reynolds  stresses  measured  by  the  CASP  will  help  model  eddy  viscosity  and  turbulence 
throughout  the  water  column.  Large  and  small  scale  bathymetry,  which  was  recorded  by  the 
CRAB  with  side  scan  sonar  and  altimeter,  and  onboard  the  sled  with  a  675  KHz  rotating  side- 
scan  sonar  and  a  2  MHz  scaiming  altimeter,  will  allow  for  a  more  comprehensive  analysis  of 
ripple  effects. 
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